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Abstract 

We consider curved space quantum corrections to the equations of motion 
of the inflaton field in the early Universe. Using the stochastic formalism in 
phase space, we demonstrate that the quantum corrected evolution of the in- 
flaton can differ dramatically from its classical evolution when the mass scales 
in the potential become large, which is naturally the case in fundamental theo- 
ries describing Planck scale physics. Using the example of the cosine potential, 
we show that the prolonged, perhaps even eternal, quantum-inflationary pe- 
riod is expected with a significant probability. This feature of Planck scale 
potentials can be dangerous, but it offers also a possibility of creating the in- 
flationary phase in string or supergravity models, where no natural realization 
of inflation has been found so far. 
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The cosmological standard model demands an unacceptable amount of fine tun- 
ing concerning the initial conditions of the Universe. A superluminal expansion of 
the scale factor, appearing during an inflationary epoch at a very early stage of the 
evolution of the Universe, was thought to circumvent these unnatural features 
However, the classical analysis of most of the proposed inflationary models shows 
that also inflation demands a fine tuning of either the initial state (new inflation) or 
the parameters of the potential (nearly every model). One of the reasons is the need 
for flat potentials coming from the fact that during inflation at least 60 e-folds must 
be attained to solve the "standard" horizon, flatness and monopole puzzles 

In the class of fundamental theories describing Planck scale physics (superstrings, 
quantum gravity) the natural scale, which sets the (nonvanishing) expectation val- 
ues and masses of various fields, is the Planck scale itself, Mp ^ 10^^ GeV. Hence, 
one has to face the problem of Planck scale potentials possibly playing, as we demon- 
strate below, a prominent role in the physics of the very early Universe. As important 
examples let us mention the so-called moduli fields arising in stringy models (they 
are allowed to have a potential when supersymmetry is broken), for instance the 
stringy dilaton which is believed to play an important role both in cosmology and in 
low-energy phenomenology and whose potential, perturbative or nonperturbative, 
remains undetermined. 

Classical analysis shows that sufficient inflation is difficult to achieve in the classi- 
cal potentials with the mass scales lying at Mp. Since one expects that the influence 
of quantum effects on the dynamics of the inflaton (by inflaton we understand any 
Planck scale fleld dominating the energy density of the Universe) and the metric 
is especially signiflcant as the energy density approaches M^, it is important to in- 
clude quantum-mechanical effects into the analysis. Attempts have been made to 
describe the fluctuating inflaton fleld during the slow roll phase by solving exphcitly 
the Fokker-Planck equation 0,10]. In the stochastic approach, which we will adopt 
here, quantum fluctuations of the inflaton are investigated in phase space. This 
method has been applied to the usual chaotic potentials in |Q. 

In this note we follow the stochastic method and reexamine the evolution of the 
inflaton in phase space for the cosine potential, which has been analysed classically 
under the name "natural inflation" 0,0. We present the results of the numerical 
analysis based on stochastic dynamics |^ with the mass parameters of the "natural" 
potential both lying at the Planck scale. We choose homogeneous universes with 
varying cosmic kinetic energies. A minimal coupling between gravity and the inflaton 
fleld is assumed. 

By means of stochastic dynamics, quantum fluctuations are given a statistical 
interpretation (quantum noise). The method takes into account the fluctuations in 
the local expansion rate due to the coupled evolution of the inflaton and the metric. 
The inflaton and its flrst time derivative (the inflaton "velocity") are decomposed 
into two parts with respect to the physical horizon [H is the Hubble param- 
eter): a long-wavelength (coarse-grained) part and a short-wavelength one. The 
quasi-classical fleld, whose evolution one follows, is formed by the coarse-grained 
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component, and the fluctuations arising in all Fourier modes are represented by the 
random force (the noise) in its equations of motion. 

The coupled evolution equations for the coarse-grained components of the infla- 
ton and the velocity v = dip/dt are given by 



v{t) = -3Hv{t) -¥'{({)) (1) 
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where the prime denotes d/dcf) . The first equation has the form of a classical stochas- 
tic equation of motion (Langevin-equation) where the random noise ri{t) includes 
the effects of the rapidly fluctuating Fourier modes. This random noise term is 
missing in the classical analysis. ri{t) is Gaussian random with zero mean and a cor- 
relation function {'r]{t)r]{t')) = 2 6{t — t'), which shows the Markov property of the 
process. In the approximation we work all spatial correlations in the random force 
are neglected to make the problem tractable (also allowing us to interpret ensemble 
averages as space averages) . The cosine potential is given in terms of the parameters 
A and / 



V{4>) = A' 



1 + cos - 



(2) 



For simplicity we use A = / = Mp, although the analysis stays valid for any set 
of Mp-magnitude scales. The coarse-grained parts of the fields, and v, are inter- 
preted as classical random variables, whose evolution is modified by the influence 
of the stochastic noise. The average over a large number of independent realiza- 
tions of the trajectories of the coarse-grained component is considered to be the 
classical homogeneous background. Treating the problem in phase space with the 
two independent phase space variables and v has the advantage that initial condi- 
tions with a dominating kinetic energy contribution, | 0^ ^ ^(0) are allowed. The 
method makes it also possible to control explicitly the fluctuations at the initial and 
final stages of inflation. There exists a well known procedure of solving a set of 
Langevin-equations by computer simulations, the molecular dynamics method 0. 
Although the accuracy of the method is not very good (it is of the order r after 
N steps of integration; r is the step length), it provides a good estimate for the 
behaviour of the mean trajectory in phase space which is finally compared with the 
classical evolution. Usually, Ito's definition of the stochastic noise is used . The 
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discrete evolution equations then take the form 



^3/2 

!n+l = Vn- SHnVn T-V^T (3) 
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where the index n characterizes the nth step of integration. The stochastic noise 
is simulated by a random number generator producing uniform deviates. They 
are turned into Gaussian deviates Wn with zero mean = and unit variance 
{wnWn') = Snn' by the Box-MuUer method 0. Simultaneously, the iterations are 
performed using a different random number at each step. Finally, the average for a 
large number of realizations is taken. 

The choice of a different interpretation of the stochastic noise is usually expected 
to have only a small influence on the final result [Q. Nevertheless, we perform the 
integration according to Stratonovich's rule to show that deviations are possible. 
In the Stratonovich case the deterministic drift term v in the first equation of the 
system (|l|) acquires an additional contribution, the noise-induced drift: 



^ noise 2 d(h 2 dd) 87r2 IQn^ ' ^ ^ 



(2) 

D)^^ is the only nonvanishing entry of the 2 x 2-difFusion matrix. Using the Stratonovich 
rule, the first equation of (|^) has to be replaced by 

^3/2 

0„+l = (t)n + VnT + L'ioLe, n^ + ^= WnV^T . (5) 



We present the results of the numerical integration of the system in the dimen- 
sionless variables 

for the cosine potential given by (1). In Fig. the evolut ion of the stochastic 
mean trajectories (according to Ito's (I) and Stratonovich's (S) interpretation of the 
stochastic noise) are compared with the classical evolution of the infiaton in phase 
space. In the example we present here, the initial conditions are chosen in such a 
way that the total energy density of the infiaton is 2 and the field is positioned 
at 0// = 37r/2. This initial value of the field is far away from the region around the 
maximum of the potential, where, by some fine tuning, one can enhance the number 
of e-folds in the classical case. No fine tuning of the initial conditions is required 
in our case. The mean trajectories shown are averages over 3 x 10^ independent 
stochastic trajectories. They share three characteristic features: (i) significant de- 
viations from the classical path occur close to the point where, they approach the 
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classical slow roll trajectory, (ii) the mean velocity tends toward zero before the os- 
cillating regime is reached, (iii) the averaged phase space evolution clearly depends 
upon the interpretation of the stochastic noise. Stratonovich's definition (S) prefers 
higher energies (this feature agrees with some earlier observations ||^). The general 
behaviour of the two mean trajectories is however similar. 
The ratio of the noise-induced drift to the classical drift, 

is also shown in Fig. Here the equation (^) has been used together with the 
equality D^^il^^ = v. The noise-induced drift acts always opposite to the classical 
one during the stage of the classical slow roll, so that the awerag'ec? evolution is being 
decelerated until the mean velocity vanishes. 

The averaged quantities are not the whole story. One has to consider the dis- 
tribution of the individual realizations around the mean. This is shown in Fig. ^ 
and Fig. |^ for the field and the corresponding velocity according to the Ito's rule. 
The numbers of integration steps are indicated. The step length was chosen to be 
T = 10~^, SO that 100 steps are equivalent to one Planck time M~^. The steps are 
also given in Fig. || for the classical trajectory. 

At the beginning, t ~ 0{M~^), when diffusion dominates over other dynamical 
effects, the distribution of the 0-values around the mean remains Gaussian in its 
shape. Gradually, the 0-dependent (nonlinear) drift coefficients cause deviations 
from the Gaussian distribution, which get significant at the time when the classical 
evolution would have already entered the epoch of coherent oscillations. It is clearly 
visible that at this stage the field has acquired values ranging over several minima 
of the potential. The maxima of the distribution are located near the minima of the 
potential. Diffusive processes are able to make the field drop into several different 
vacua. 

The f-distribution is also purely Gaussian at the beginning although this is not 
shown in Fig. 0. The nonlinearity of the potential gives rise to significant deviations 
from the Gaussian shape already long before the initial velocity has decayed. It 
is remarkable that the distribution of this variable differs from Gaussian statistics 
considerably earlier than that of the 0-variable. Gradually, it develops a tail toward 
positive velocities. At this stage the mean velocity constantly decreases, whereas 
the tail forms a second maximum, so that the distribution becomes symmetric with 
respect to the mean. This is not unexpected since the quantum "kicks" appear 
with similar probability in opposite directions (in the 0-space), so that upward and 
downward motions are generated in equal abundances. One should note that the 
total width of the ^-distribution is much smaller than that of the 0-distribution. 

The final stage of the stochastic evolution is characterized by a zero mean veloc- 
ity. This is due to the symmetry of the w-distribution. The averaged drift is zero, 
so that the averaged field seems to stick to one value. However, the most probable 
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values of the velocity are by no means zero (thus indicating the dominance of the 
fluctuations in the inflaton dynamics); positive velocities are then equally proba- 
ble as negative velocities, so that the noise-induced movement uphill has nearly 
the same probability as the downward motion. We have found differences between 
the time evolution of the inflaton for Ito's and Stratonovich's interpretations of the 
stochastic noise. But both cases share the same basic feature: since the noise- 
induced drift always acts against the classical drift, the time evolution of the mean 
quasi-classical field (compared with the classical case) can be delayed by quantum 
fluctuations. The stochastic analysis shows that the duration of the quasi-de Sitter 
phase can be prolonged to reach a large number of e-folds, much larger than 60. 
The a priori probability for this to happen is estimated to be of the order of one 
(see discussion below). 

The behaviour described so far is generic for any Planck scale field dominating 
the energy density of the Universe. The discussion demonstrates clearly that the 
classical description of the cosmological evolution driven by such a field is highly 
inadequate and that a prolonged quasi-de Sitter phase is naturally generated in the 
quantum description. 

Let us try to understand whether this quantum inflation can provide us with a 
model for successful inflation, solving cosmological puzzles. The basic requirement 
for the scenario under consideration to serve as a successful inflationary model is 
that the number of e-folds produced during the superluminal stages of the evolution 
of the Universe must be large enough, larger than 60 at least. Of course, there are 
model universes in the ensemble where the inflaton rather rapidly approaches values 
near the minimum of the potential and then stays there. In this case the generated 
number of e-folds is insufficient. 

However, there is a large enough fraction (of order unity) of all observed realiza- 
tions of the scenario where the fluctuations act in such a way that the field varies 
(over macroscopic periods of time) much more slowly than required by the classical 
evolution. In this case, a sufficiently large number of e-folds could be generated 
to solve the cosmological puzzles. Fig. ^ shows the development of the number 
of e-folds produced during the evolution of the inflaton for one typical stochastic 
path compared with the classical result. Apparently, the quantum dominated mo- 
tion easily provides more than 60 e-folds. This means that an inflationary epoch is 
generated by quantum fluctuations for a generic inflaton potential with parameters 
fixed at the Planck scale. Our results indicate that even the a priori probability 
of quantum inflation is large, not to mention the a posteriori probability , which 
is exponentially close to one. We stress again, that the prolonged quasi-de Sitter 
stage is a generic feature of our scenario and does not require a fine tuning of the 
initial conditions. 

The end of the inflationary period and the reheating of the Universe in this 
scenario have to be examined more carefully. Generally, the strength of the fluc- 
tuations increases with increasing values of the potential (as the Hubble parameter 
increases), and the field can easily be driven uphill thus leading to a prolongation 
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of the quasi-de Sitter phase. Hence, the possibihty of "eternal" inflation is in- 
corporated in this model. A natural exit is always possible with the assumption of 
the instability of the infiaton. By decaying into light particles a sufficient reheating 
temperature necessary for baryogenesis could be provided and inflation be stopped. 
A specific feature of this violently fiuctuating system is the non-zero probability 
of large fiuctuations throwing the field out of the infiationary regime thus ending 
abruptly the quasi-de Sitter phase. 

One possible problem of quantum Planck scale infiation is the excessive produc- 
tion of energy density inhomogeneities. Since the relative magnitude of a physical 
length scale to the horizon, X/H~^ oc i?, is a fiuctuating quantity, any scale A crosses 
the horizon not only once, but frequently within the infiationary period. For this 
reason, the resulting amplitude of generated density perturbations need to be cal- 
culated carefully, as it would serve as a very sensitive probe of the present scenario. 
One obvious way out of the density fiuctuations problem is to assume that the final 
stage of infiation, when density perturbations at observable scales are produced, is 
controlled by some other, lower-scale, field (like in multiple-infiaton models [0), 
which can easily be the case in models, like string inspired ones, where the number 
of relevant degrees of freedom is rather large. In such a case, quantum infiation is 
the natural trigger for the subsequent infiationary epoch and can easily provide a 
quasi-de Sitter stage of realistically long duration. 

On the basis of our experience with a wide class of potentials (we have found 
similar effects for the quadratic and quartic potentials of chaotic inflation) we want 
to conclude that the presence of quantum inflation seems to be a generic feature of 
fleld theoretical models addressing Planck scale physics and it may be interesting 
and fruitful to explore this phenomenon in more detail. 
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Figure 1: Stochastic mean trajectories (averaged over 3 x 10"^ realizations) in the 
potential V{(f)) = [l + cos(0//)] according to Ito's (I) and Stratonovich's (S) 
interpretation of the stochastic noise. The parameters A and / are set to the Planck 
scale. The classical trajectory is shown by the dotted curve. It meets the slow roll 
curve (long dashed) after about 30 steps of integration. The initial condition is 
chosen to lie on the curve of constant energy density p = 2 A^ (dot-dashed hues) 
with a field value of 0// = 37r/2 and positive velocity. 3 x 10^ integration steps were 
performed with a step size of r = 10~^. Several hnes of constant ratio TZ (dashed 
lines) are shown. The absolute value of 71 increases, if one moves towards ij = 0. 
The values of 71 and the number of integration steps along the classical trajectory 
are indicated. 
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Figure 2: Distribution of the 0-values of 10^ realizations (normalized) around the 
mean (x=0) after 100, 1000 and 5000 steps of integration. The influence of the 
potential is visible. 
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Figure 3: Distribution of the t;-values of 10^ realizations (normalized) around the 
mean (x=0) after 100, 300 and 5000 steps of integration. 
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Figure 4: Total number of e-folds produced during the evolution of the inflaton. The 
initial condition and the parameters of the potential are chosen as in Fig. 1. The 
classical analysis (C) clearly shows the end of inflation after the field has reached the 
minimum of the potential. Only 25 e-folds can be produced before the entry into 
the hot Friedmann universe. The time behaviour of a typical stochastic solution (S) 
reveals a continious superluminal expansion of the universe; the demanded number 
of 60 e-folds is easily exceeded. 
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